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P_i ' We present an approach to the canonical quantization of systems with equations of motion that are 

historically called non-Lagrangian equations. Our viewpoint of this problem is the following: despite the 
^ fact that a set of differential equations cannot be directly identified with a set of Euler-Lagrange equations, 

one can reformulate such a set in an equivalent first-order form which can always be treated as the Euler- 
Lagrange equations of a certain action. We construct such an action explicitly. It turns out that in the 
^ ■ general case the hamiltonization and canonical quantization of such an action are non-trivial problems, since 

' the theory involves time-dependent constraints. We adopt the general approach of hamiltonization and 

canonical quantization for such theories (Gitman, Tyutin, 1990) to the case under consideration. There 
exists an ambiguity (not reduced to a total time derivative) in associating a Lagrange function with a given 
set of equations. We present a complete description of this ambiguity. The proposed scheme is applied to the 
quantization of a general quadratic theory. In addition, we consider the quantization of a damped oscillator 
and of a radiating point-like charge. 



1 Introduction 



It is well-known that some physical systems like dissipative systems [I], Dirac monopole [2], etc. are usually 
described in terms of second-order equations of motion which cannot be directly identified with Euler-Lagrange 
equations for an action principle. Following traditional terminology, we call such equations of motion non- 
Lagrangian equations in what follows. Sometimes (but not always) non-Lagrangian equations can be reduced 
to Euler-Lagrange equations by multiplying by the so-called integrating multiplier [3]- [5]. The existence of an 
action principle for a given physical system, or what is the same, the existence of a Lagrange function for such a 
system, allows one to proceed with canonical quantization schemes. This, in particular, stresses the importance 
of formulating action principle for any physical system. 

In the present work we discuss an approach to constructing quantum theories that in the classical limit 
reproduce non-Lagrangian equations of motion for mean values. In fact, we consider a canonical quantization 
of Lagrangian theories with time-dependent constraints that are related to the non-Lagrangian systems. To 
this end, on the classical level, we reduce non-Lagrangian equations of motion to an equivalent set of first-order 
differential equations. For such equations, one can always construct an action principle, the corresponding 
consideration is represented in section 2 and, partially, is based on results of works [8]-|10|. The hamiltonization 
of the constructed Lagrangian theory leads to a Hamiltonian theory with time-dependent constraints as it is 
demonstrated in section 3. Thus, we show that systems traditionally called non-Lagrangian ones are, in fact, 
equivalent to some first-order Lagrangian systems, however, with time-dependent constraints in Hamiltonian 
formulation. The canonical quantization of the latter theory is not a trivial problem (it follows to general 
consideration of [14) ) and is represented in section 4. It is known that on the classical level, there exists an 
ambiguity in constructing Lagrange function (which is not reduced to a total time derivative) for a given set of 
equations [8]- [13]. We describe completely such an ambiguity for the case under consideration. We apply the 
general approach to formulate the canonical quantization in case of theories with arbitrary linear inhomogenous 
equations of motion (general quadratic theories), see section 5. Then we consider the canonical quantization of 
a damped harmonic oscillator (sec. 6) and a radiating point-like charge (sec. 7). 
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2 Action principle for non-Lagrangian systems 

Let a system with n degrees of freedom be described by a set of n non-Lagrangian second-order differential 
equations of motion. To construct an action principle, we replace these equations (which is always possible) by 
an equivalent set of 2n first-order differential equations, solvable with respect to time derivatives. Suppose such 
a set has a form 

x a = f a (t,x), o=l,..,2n, (f) 

where f a (t,x) are some functions of the indicated arguments and by dots above we denote time derivatives of 
coordinates. Since these equations are first-order, action S[x] that yields (QJ as Euler-Lagrange equations, must 
be linear in the first time derivative x a . Its general form is 

S[x]= JdtL, L = J a x a -H 1 (2) 

where J a = J a (t,x) and H = H(t,x) are some functions of the indicated arguments. The Euler-Lagrange 
equations corresponding to ((2]) are 

«!^-^!§ = 0=^ -d a H - d t J a + (d a J p - d p J a ) x^O, (3) 
where the notation are used 

o 9 a - 9 
Denoting the combination (d a Jp — dpJ a ) by £l a p , 

= da J/3 ~ dp J a = Cl a/3 (t,x) = -£l/3 a (t,x) , (4) 

we rewrite (jSJ as follows: 

n aP ±v = d a H + d t J a . (5) 

Equations (jSJ or ||SJ| can be identified with (TIJ, provided 

dot n 0/s ^ a , (6) 

Va/sf ~ d t J a = d a H . (7) 

The functions J a and H can be found from conditions if the matrix Q a p is given. Assuming that J a 

and H are smooth functions the consistency condition for equations ([7]) imply 

d («a 7 / 7 ) ~ d a (SlfrP) + d t n aP = o => d t n aP + £ f tt afi = , (8) 

where £ f£i a /3 is the Lie derivative of £l a /3 along the vector field f 1 . In addition, one can verify that the matrix 
O obeys the Jacobi identity (Q a /3 is a symplectic matrix) 

d a flp-y + dpVt ia + d 7 n a f} = . (9) 

Now we are going to analyze these equations. It is known that the general solution Cl a p of equation ([8]) 
can be constructed with the help of a solution of the Cauchy problem for equations ([!}. Suppose that such a 
solution is known, 

x a =cp a (t, X{0) ), xf Q) =cp a (0,x (0) ) (10) 

be a solution of equations (TT]) for any x(o) = (^(o)! ' ano - X a {ti x ) be the inverse function with respect to 
P a (*>ff(o)), i-e., 

x a = v a (t,x (0) )^x? 0) = X a (t,x), x a = ^ a (t, X a ), 0aX 7 |t=o = ^. (11) 

Then 

n a p(t,x) = d aX ' 1 ci i °l(x) dp X 5 , (12) 



where the matrix Q^l is the initial condition for f2 Q ^ 



ftap(t,x)\t=o = ^-a}( x ) 
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It follows from (J9j) at t — that the matrix Cl^l(x) obeys the Jacoby identity such that the general structure 
of this matrix is (we do not consider global problems which arise from a nontrivial topology of the x"-space) 

n i°3 = daj/j - d j a , (13) 

where j a (x) are some arbitrary functions. Then equation (| 1 2[) implies 

= ~ 9ptp a > ^a(*) 2; ) = J> (x(*i »)) ^aX^ '(*) &) ■ ( 14 ) 

On the other hand, relation (j4]) must hold, 

- 9pip a = d a J/3 - dpJ a , 

which implies that 

J a (t, x) = tp a + d a f = jp (x(t, xj) d a X (t, x) + d a tp(t, x) , (15) 

where tp(t,x) is an arbitrary function. One can represent another form for J a (t,x), in which the ambiguity 
related to the arbitrary functions jp (x) is incorporated in the matrix £1^1 • To this end, we remind that the 
general solution for J a (t,x) of the equation |4|) provided that Q a p is a given antisymmetric matrix that obeys 
the Jacobi identity, is given by 



J a (t,x) — I xPfHp a (t, sx) sds + d a (p(t,x) , (16) 
Jo 

where ip(x) is an arbitrary function. Substituting (|12j) into l|16p . we obtain 

J a (t, y) = f y fi kx 7 V$ (x) dpx 5 } sds + dMt, y) ■ (17) 



Equations (|15p or (??) desribe all the ambiguity (arbitrary functions jp (x) and ip(t,x), or arbitrary symplectic 

S<5 



matrix and arbitrary function (p(t, x)) in constructing the term J a (t, x) of the Lagrange function 



One can also see that choosing the matrix rt^(x) to be nonsingular, we guarantee the nonsingularity 
(condition ©) for the matrix £l a p(t,x) since components of the latter are given by a change of variables (TT2)) . 

To restore the term H in the Lagrange function ([2]), we need to solve the equation with respect to H. 
To this end, we remind that the general solution of the equation dif = gi, provided a vector is a gradient, is 
given by 

f(x) = ds x z gi(sx) +c, 
Ja 

where c is a constant. Taking the above into account, we obtain for H the following representation: 

H(t,x)= [ dsx? [tlp a (t,sx)f a (t,sx) - d t Jp(t,sx)]+c(t) , (18) 
Ja 



where c(t) is an arbitrary function of time, and flp a and Jp are given by (| 12(1 and (|17p respectively. All 

)(0) 
Si ' 



the arbitrarness in constructing H is thus due to arbitrary symplectic matrix 0^$ , arbitrary functions (p(t,x) 



entering into VLp a and Jp and due c(t). 

We see that there exist a family of actions ((2|) which lead to the same equations of motion (jTJ) . It is easy 
to see that actions with the same but different functions ip(t, x) and c(t) differs by a total time derivative 
(we call such a difference trivial). A difference in Lagrange functions related to different choice of symplectic 
matrices 0,^1 is not trivial. The coresponding Lagrangians are known as s-equivalent Lagrangians. In spite of 
the fact that actions with nontrivial difference lead to the same equations of motion, they lead in general to 
different Hamiltonian formulations and to different quantum theories in course of the quantization. However, 
any quantum theory that is obtained by the developed below quantization procedure obeys the correspondence 
principle, i.e., in the classical limit, equations of motion for mean values coincide with ([1]). Equations of 
motion |T]) do no contain any additional information that can be used in choosing a "right" quantum theory. 
Only physical considerations or a comparison with experiment may be used for this aim. Below, we are going 
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to consider hamiltonization and subsequent quantization of the action ^ with the following choice of the 
symplectic matrix f2^ 

< 19 > 

where I is an n x n unit matrix, and denotes an n x n zero matrix. This choice of the action leads to the 
canonical commutation relations for original variables on the quantum level. Hamiltonization and quantization 
of the action @ with different choices of the symplectic matrix ilty can be fulfilled in the same manner, but 
technically look more clumsily. 

Note that (JTUJ) implies that there exist only two possibilities for the matrix £1 in (fT2"| . Namely, it is either 
a canonical symplectic matrix, which is possible only if the initial equations ([1]) are canonical Hamiltonian 
equations (or, equivalently, Lagrangian equations of motion for the first-order action), or it must depend on 
time, which is the case of non-Lagrangian equations. 

The first-order action ([2]) can be regarded as a Lagrangian action, or as a Hamiltonian action with a 
noncanonical Poisson bracket. An equivalent second-order Lagrangian formulation is always possible; however, 
it may include additional variables [15] . 

One ought to say that it is always possible to construct a Lagrangian action for non-Lagrangian second-order 
equations in an extended configuration space following a simple idea first proposed by Bateman 22J. Such a 
Lagrangian has a form of a sum of initial equations of motion being multiplied by the corresponding Lagrangian 
multipliers, new variables. Euler-Lagrange equations for such an action contain besides the initial equations 
some new equations of motion for the Lagrange multipliers. In such an approach one has to think how to 
interprete the new variables already on the classical level. Additional difficulties (indefinite metric) can appear 
in course of the quantization. 

As an example, we consider a theory with equations of motion of the forrrjj 

x = A(t)x + j(t) . (20) 

We call such a theory the general quadratic theory. Let us apply the above consideration to construct the action 
principle for such a theory. 

Solution of the Cauchy problem for the equations (|20|) reads 

x(t)=T{t)x {0) + 7 (t), (21) 

where the matrix T(t) is the fundamental solution of ([20]) . i.e., 

f = AT , r(0) = 1 , (22) 

and 7(t) is a partial solution of (j20|) . Then following (fT2"f . we construct the matrix f2, 

q = a t q^a, A = r _1 . (23) 

and find the functions J and H according to (fT"7)) and (fT8|) , 

J = -xVL , H = -xBx - Cx , (24) 
2 2 

where 

B = - (OA- A T fl) , C = Qj. (25) 
Thus, the action functional for the general quadratic theory is 

S[x] = \f dt ( xn± ~ xBx ~ 2Cx ) ■ ( 26 ) 

Another approach to constructing the action functional for the general quadratic theory was proposed in [T] . 
Note that Darboux coordinates xq can be written via a matrix A as follows: 

x -> x = R^i^A^x. (27) 



l Here we use matrix notation, x = (x a ) , A(t) = \ A{t)^\ , j(t) = , o,/3 = 1, ..,2n. 
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Here, R(t) is an arbitrary matrix of a linear (generally time-dependent) canonical transformation: 

R T (t)Q^R(t) = n<°> . 
In terms of the coordinates xq , action (J^U) takes the form 

S[x] = \f dt (x n {Q] x + XoR T n^Rxo - 2CTRx^j . (28) 

The Darboux coordinates (|2"T)) can be divided into coordinates and corresponding momenta. The Euler- 
Lagrange equations for action J28|) have the form of canonical Hamilton equations with the Hamiltonian 

H = - ^x R T n {0) Rx Q + CTRx . (29) 

Note that the choice R = const yields a trivial Hamiltonian, which is consistent with the fact that in this 
case Xq are the initial data without dynamics. 

3 Hamiltonian formulation 

We are now going to consider action ([2]) as a Lagrangian action with the Lagrange function 

L = J a (t, x)x a -H (t, x) (30) 

and construct a corresponding Hamiltonian formulation. To this end, we follow the general! scheme of |14) . We 
first construct the action S v [x,ir,v], which, in this case, has the form 

S v [x,tt,v}= J [J a (t,x)v a - H(t,x)+7r a (x a -v a )]dt, (31) 

and depends on the momenta 7r Q conjugate to the coordinates x a , as well as on the velocities v a . The equations 

5S V 

— = ® a (t,x,TT)=n a -J a (t,x)=0 (32) 
ov a 

do not allow one to express the velocities via x and 7r, which implies the appearance of primary constraints 
$ Q (t, x, 7r), and the velocities v a become Lagrangian multipliers to these constraints, so that action (f3Tj) becomes 
a Hamiltonian action of a theory with the primary constraints (|32[) , 



Sh = J dt{w a x a - , = H(t, x) + A Q $ Q (t, x, tt) , (33) 

with the equations of motion 

f,= {r],H^] , $ = 0, (34) 

where rj = (x,tt). 

The primary constraints are second-class ones. Indeed, we have, in virtue of (O, 

$^1 = n aP (t, x) => dct{$ Q , ^1^0. (35) 

Thus, secondary constraints do not appear, and all A-s arc determined from the consistency conditions for the 
primary constraints: 

$ Q = d t <S> a + {$ Q , = => -d t J a - d a H + A^{$ Q , <S> fj } = 

\ =u 0a {d t J a + d a H) , = fi^. (36) 



2 Note that some of J a can be equal to zero, for instance, if one deals with a canonical Hamiltonian action. In this case, one 
obtains the constraints <J> a = n a = 0. Another way to examine this case is to use the method of hamiltonization for theories with 
degenerate coordinates 1161 . 
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Using the Lagrange multipliers (|36|) in equations (|34|) . we can write these equations in the form 

r) = {q, H} Dm + {77, $> a }u a0 d t J? , $ = , (37) 

where {• • • , • • • ar e the Dirac brackets with respect to the second-class constraints $. For the canonical 

variables the Dirac brackets are 

{x a ,x l3 } Dm 
{ttqjTt p}d(<s>) 

{x a ,TTf3} D ^) 

Formally introducing a momentum e conjugate to the time t, and defining the Poisson brackets in an extended 
space of the canonical variables (x, tt; t, e) = (7/; t 7 e), see [H], we can rewrite ([3"7|) as follows: 

f] = {r h H + e} Dm , $ - . (39) 

Equations (|39p present a Hamiltonian formulation of non-Lagrangian systems with first-order equations of 
motion |T]). We note that the Hamiltonian constraints in this formulation are second-class ones and depend 
on time explicitly. The canonical quantization of theories with time-dependent second-class constraints can be 
carried out along the lines of [14] . Below, we present the details of such a quantization, and then adopt it to 
the system under consideration. 



a/3 



d a J p u) pl df3J~ l , 

5% + UJ^dgJ, 



(38) 



4 Canonical quantization 

For a Hamiltonian theory with time-dependent second-class constraints, the quantization procedure in the 
"Schrodinger" picture is realized as follows. The phase-space variables 77 of a theory with time-dependent 
second-class constraints $i(rj,t) are assigned operators fj (t) subject to the equal-time commutation relations 
and the constraints equations 

{v A (t),V B (t)}=i{ V A , V B } DW \ ri=fl , $i(?i(t),t) =0. (40) 

Their time evolution is postulated as (we neglect the problem of operator ordering [17] ) 

d 



dt 



V (*) = {V; \n=n = ~{»7. ^H $ ' ®}iv d t®i> \v=v ■ ( 41 ) 



To each physical quantity F, given in the Hamiltonian formulation by a function F{t, 77), we assign a "Schrodinger" 
operator F (t) , by the rule F (t) = F(t, 7} (t)). For arbitrary "Schrodinger" operators F (t) and G (t) , the relation 

[F(t),G(t)]=i{F,G} DW \ v=fl (42) 

holds as a consequence of (I4U1) . The quantum states of the system are described by vectors $ of a Hilbert space 
with a scalar product ("J, V E ,/ ). Their time evolution is determined by the Schrodinger equation 

= (43) 

where the quantum Hamiltonian H is constructed according to the classical function H (t, rf) as H (t) = 
H(t, 77 (t)). The mean values (F) t of a physical quantity F are determined as the mean values of a corresponding 
"Schrodinger" operator F (t) = F(t, fj (t)) with respect to state vectors ^ (t) , 

(F) t = (*(t),F(t)V(t)) . (44) 

Provided that H is a self-adjoint operator, the time evolution of state vectors \1/ (t) is unitary, 

* (£) = u (t) * (0) , U + (t) = U- 1 (t) , (45) 
where U (t) is an evolution operator. 
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In the Heisenberg picture, where state vectors are "frozen" and the time evolution is governed by the 
Heisenberg operators f] (t) = U^ 1 (t) f] (t) U (t), one can see [H] that 

—fj = {r), H(t, rj) + e} D( ^ \ v=f , , 

[f, A (i) , ff (t)} = i{r, A , V B } Dm \ v=f> , Hf, (t),t)=0, (46) 
while for Heisenberg operators F (t) = U^ 1 (t) F (t) U (t) = F{t, fj (t)), we have 

j f F (t) = {F(t, r)),H(t, n) + e} DW \ v =t , (47) 

or 

j t F (t) = -i [F (t) , H (t)] + {F(t, V ), e} Dm \ n = n . (48) 

The mean values (F) t in the Heisenber picture in according to (|4"4"]) and ([43]) are determined as 

(*%=(*(0),F(t)¥(0)) . (49) 

The above quantization provides the fulfilment of the correspondence principle because quantum equations 
has the same form as the classical one 

Note that the time-dependence of the Heisenberg operators in the theories under consideration is not unitary 
in the general case. In other words, there exists no such ( "Hamiltonian" ) operator whose commutator with a 
physical quantity can produce its total time derivative. This is explained by the existence of two factors which 
determine the time evolution of a Heisenberg operator. The first one is the unitary evolution of a state vector 
in the "Schrodinger" picture, while the second one is the time variation of a "Schrodinger" operators f], which 
in general has a non- unitary character. The existence of these two factors is related to the division of the 
right-hand side of (|4"8|) into two summands. Physically, this is explained by the fact that dynamics develops on 
a surface which changes with time - in the general case, in a nonunitary way. 

Below, we apply the above quantization scheme to the system under consideration. Taking into account the 
Dirac brackets (|55|) , we can write the equal-time commutation relations (|4U|) for phase-space operators as 

[TTa,^] = i daJpW 91 dp,J-f\ x=s . , (50) 

[£",*{,] =i6% + iu aif df l J J \ x=ji . 

In this case, the classical Hamiltonian H does not depend on the momenta n a , and therefore in order to 
determine the quantum Hamiltonian H, we need to know only the time dependence of the operators x a . From 
(|4l)l it follows that 

±x a = u a P(t,x)d t J f} (t,x)\ x=i . (51) 

5 Quantization of general quadratic theory 

The quantum-mechanical description of quadratic systems is a widely discussed physical problem which has a 
number of important applications (see, e.g., [H]-[21] and references therein). Almost all of these works deal 
with the case of "Hamiltonian" quadratic systems, i.e., systems described by canonical Hamiltonian equations 
of motion. On the other hand, we consider a general quadratic system, i.e., a system described by arbitrary 
linear inhomogeneous equations of motion ([20| . In this case conditions (f50l) , (f5T]l become 



[x a ,x p ] =itJ aP (t), (52) 

j t x a = ~ l -^m p ,{t)x\ (53) 

The time-dependence of the operators x can be easily found: 

x a (t) = 3> a P(t)x$. (54) 
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Here, the matrix $ obeys the equation 



$ = -^cuti®, $(0) =E, (55) 



and the operators Xq obey the following commutation relations: 

Xq , Xq 



(nS3)"' = <( _°, J )• <») 



see (|19l) . In what follows, it is useful to divide the operators Xq into the operators of coordinates proper and 
corresponding momenta Xq — (q l ,pi), a — l,..,2n, i = l,..,n. The operators q and p obey the canonical 
commutation relations 

[q\pj] =i6), [q\qi] = = 0. (57) 

The quantum Hamiltonian in eq. (|43[) takes the form 

H = ^x $ T B$x - C$io , (58) 

where the matrix B is determined by (125j) . 

The above quantization is equivalent to quantization in Darboux coordinates, and the transformation x — > 
$(t)xo provides, by itself, a passage to the Darboux coordinates xq, because (|55[) implies 

= fl . (59) 

Namely, in the coordinates xo the Poisson bracket is canonical. Therefore, $ = T(t)R(t), where T(t) is a 
fundamental solution of system (|2U)) . However, in contrast to the classical theory, now the matrix R(t) is fixed, 
it must obey the conditions 

R = fl T t BTR, R(0) = E. (60) 
Thus, using (|29|) one can also rewrite the Hamiltonian in (f58|) as follows: 

H = --xoR T Sl^Rx + CTRx ■ (61) 

It is remarkable that if the matrix A that determines the set of equations ([2T)]) is constant, the matrix that 
determines the quadratic part of the Hamiltonian in (|58p is constant as well, and equals to 



$ T B$ = 5(0) = ^(o'°)A-A T l]( >). (62) 

This fact is easy to observe because the time derivative of this matrix, in view of (f55")) . (f2"5|) and ([23)) , is equal 
to zero: 

— ($ T S$) = 0. 

Thus, in this case, as distinct from the general case, the matrix $ can be determined from the set of algebraic 
equations (f5"9"]) and (p2^) . 

Note that if we start from a canonical Hamiltonian system the above quantization coincides with the usual 
canonical quantization, because in this case equation (|53[) becomes dx/dt = 0, i.e., x(t) = Xq. 
In the Heisenberg picture, equations (|46[) for the operators x take the form 

±x = A(t)x + j(t), (63) 
[x a ,x fs ]=iuj afs {t) . (64) 

Equations (|63p coincide (the correspondence principle) with the classical equations of motion (|20p : however, the 
commutation relations (|64|) differ from the canonical ones. So, evolution of operators x can be written as 

x (t) = T(t)x + 7 (t) , (65) 
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where operators Xq as well as xq obey the canonical commutation relations 

= '(««) "' = '(-! »)■ m 

Thus, mean values (F)t of a physical quantity F according to (|49j) are determined as mean values of the 
corresponding operator F (t) = F(t, x (i)) with respect to initial states vectors (0) , i.e. 

(F) t = (* (0) , F {t, T(t)x + 7 (t)) * (0)) . (67) 

We see that the quantum evolution of physical quantities in general quadratic systems is completely determined 
by the classical one. 



-a ~P 

Xq , Xq 



6 Quantization of damped harmonic oscillator 

The above formulated quantization of non-Lagrangian theories and, in particular, of general quadratic theories 
can be immediately applied to quantizing a damped harmonic oscillator. The latter problem attracts attention 
for already more then 50 years , there exist different approaches to its solution, no one of them seems to be a 
final version which does not contain weak points, see e.g. [23] -02], PQ. 

The classical equation of motion for a damped harmonic oscillator is non-Lagrangian, it has the form 

r + 2ar + uj 2 r = 0, (68) 

where w is the angular frequency and a > is a friction coefficient. Introducing an auxiliary variable y = f, we 
reduce (|68[) to the following equivalent pair of first-order equations: 



r = y , y = —u 2 r — 2ay . (69) 
Following the way proposed in sec. 2, we construct an action S that implies ([6"9"]> as Euler-Lagrange equations, 

S = - J dt[yr-ry- (y 2 + 2ary+uj 2 r 2 )] e 2at . (70) 

Note that equation ([68| can be represented as 

-(e 2at r)+e 2at u; 2 r = 0, 
at 

i.e., as a Lagrangian equation of motion with time-dependent mass and frequency. In this case, the mass e 2at is 
nothing else but an integrating multiplier for equation (|68[) ; however, as was already mentioned, an integrating 
multiplier does not always exist [3J [TO] . 

Then we proceed with the canonical quantization described in the previous section. Equal-time commuta- 
tions relations (|52|) and equations (|53[) determining time evolution of "Schrodinger" operators r and y are 

[r,y] = ie- 2a \ [f, f] = [y, y] = Q, (71) 

P = " af > iy = - ai) - (72) 

A solution of these equations has the form 

r = e- at q, y = e~ at p , (73) 
where operators q and p obey canonical commutations relations 

[q,P] = i, = \P,P] =0. 
According to (|58|l . the corresponding quantum Hamiltonian reads 

H = - [p 2 + a(q P + pq)+uj 2 q 2 ] . 
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P 2 + (to 2 - a 2 ) Q 2 (74) 



It can be modified to the form 

H=l 
2 

by the help of the canonical transformation (p, q) — > (P, Q), where P = p + aq, and Q = q. The corresponding 
generating function is = qP ~ aq 2 /2. 

As usual we define the classical energy of the system by 

E= l -(r 2 +oj 2 r 2 )=\(y 2 +uj 2 r 2 ) . 

One can easy see the energy depends of time as follows: E — Eoe~ 2at . Using (|73p. we obtain an expression for 
the operator E that corresponds to the classical quantity E, 



E = -e- 2at 
2 



P z -a{PQ + QP) + (Lu 2 + a 2 )Q 2 . (75) 

Let us consider the underdumped case, a < u>. Then (|74|) is a Hamiltonian of a harmonic oscillator with an 
angular frequency Co — y/u; 2 — a 2 . Stationary states of the corresponding Schrodinger equation have the form 

*„ = e- lE ^ n (Q) , E n = Q (n+ , n = 0, 1, 

tf.(Q) = ^= (|) V4 e-^/ 2 F„ (QV5) , (76) 
The mean values of energy (of the operator (|7Sf ) in such states can be easily calculated, 

< £ >" = (" + 5)^ 2 °'' (77) 

At each fixed time instant, the energy spectrum is discrete, however, it decreses with time exactly as in classical 
theory. The same conclusion was derived in |36j-[38j where a second-order action obtained by integrating- 
multiplier method was taken as a starting point for a quantization. A quantization of the damped oscillator 
following Bateman (see above) meets serious difficulties such as indefinite metric etc., [381139], 

Overdumped cases, when a > u> correspond to an aperiodic motion in classical theory [43] . Its quantum 
interpretation is not clear due to continuous character of the Hamiltonian spectrum. 

A nontrivial generalization of equation (|68p could be an n-dimensional damped oscillatoJl 

r + 2ar + ur = , (78) 

with the matrixes a and us being constant and symmetric. Introducing auxiliary variables, y = r, y = [y i j , we 
reduce (|78|) to the following set of first-order equations: 

x = Ax, [x= (x a ) = (rW)), (79) 
I 



.4 



— io —2a 



This is a set of linear equations with constant coefficients. In this case, the corresponding quantum Hamiltonian 
can be constructed, according to (f6"2"]) . as follows: 

H=h o B(0)x o , (80) 

where 

/ u> a 



B(0) = 



a I 



and the operators Xq can be divided into the operators of the coordinates proper q and those of the corresponding 
momenta p, with the canonical commutation relations (|57|) . Further solution of the quantum problem with 
quadratic Hamiltonian ([50)1 can follow, for example, jT5J ?]. 



3 Here, we use matrix notation, 



r = (r 1 ) , a = \ a % A , ui = (u^j , i,j = 1, 
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7 Quantization of radiating point- like charge 



Equations of motion for a nonrclativistic particle moving in electric E and magnetic H fields with account taken 
of the back reaction of the radiation emitted by the particle have the form [H] 

mr = F + f , (r = (x,y,z)), 
e 2e 2 

F = eE+-[f x H], f =— f. (81) 
c 3c 3 

Here F is the Lorentz force, f is the force of the back reaction of the radiation, e is the charge of the particle, 
and c is the light velocity. Derivatives with respect to time are denoted by dots above. 

These equations are of third order, therefore a trajectory of a charged particle cannot by uniquely specified 
only by initial position and velocity of the particle. It was also pointed out that together with physically 
meaningful solutions, equations (|8ip have a set of nonphysical solution [44] . However, in the case when the back 
reaction force f is small when compared to the Lorentz one F, 

|f| « |F| , (82) 

these equations can be reduced to the second order equations by means of a reduction of order procedure. Then 
the above mentioned problem with nonphysical solutions does not appear. In the reduction procedure, equations 
(f5Tj) are replaced by second-order equations f = g(r, r, e) such that all the solutions of the latter equations would 
be solutions of (|8ip . The last requirement implies a partial differential equation on the function g(r, r, e) having 
a unique solution with the natural condition g(r, r,0) = 0, see e.g. [4"4l 110) . 

Consider, for example, a particular case E = 0, H = (0,0, H = const). In such a case, the reduced 
second-order equations have the form [lOj 

x = —ax — f3y , y = (3x — ay , z = , 

V6V3 + V9 + 64eSiP 6 2 effy^ 

a = —5 Rj — e H , p = = « eH . (83) 

8e 3 ^3 + V9 + 6Ae 6 H 2 

Here we have set m = c = 1 for simplicity. Since the evolution along the z-axis represents the free motion 
and decouples from the dynamics in the xy-plane, we restrict our consideration to the first two equations. At 
a = 0, equations are Lorentz equations with an "effective" magnetic field (3 = (0,0, P/e). In this case, the 
trajectories are concentric circles. If a ^ 0, the particle spirals at the origin of ccy-plane. So, it is natural to 
treat a as a friction coefficient. 

In order to construct an action functional for the non-Lagrangian second-order equations (|83p , we introduce 
new variables as follows: 

. , j8 • P 

p^ x + -^y, 9 = y- 

In the new variables, we have a set of first-order equations, 

P . ^P 
x = p- -y, y = q+-x, 

P P 2 ( P \ . P P 2 f.P \ (RA , 
V = —x- a\p- -y\ , q= -p- —y~ a\q+ -x\ . (84) 

According to general formulas (|24| . ([25]) . and ([25| . we construct the following action for the set ([84]) : 

1 

2{a~+ f) 

+ c{pq - qp) + d(xy - yx) - e(p 2 + q 2 ) - f{x 2 + y 2 ) - g{px + qy) - j{qx - py)] dt , 
where time-dependent functions a, b, c, d, e, /, g, and j are 

a = a 2 cos(/3i) + /3 2 cosh (-at) , b = a 2 sin(/3i) + a(3e~ at - a/3 cos (/ft), 



S = ^7~i ^27 / e " ^(P* - x P + <iy- yo) + P(qx -xq + yp- py) (85) 



- at f3 - 2a sin(pY), d = -y sinh (-at) - ^a(3sin((3t) + a 2 (3 [cos(/3i) - < " ' 



+ a 2 cos(/3t) - sin(/3t), g = a ,9 cos(/3t) + a J cos(/3t), 
/ = | {/3 3 e at + /3a 2 cos(/3t) + a[/3 2 + 2a 2 ] sin(/3t)} , j = a 3 sin(/3t) + e at (3 3 + a 2 /3 cos(^). 



11 



In the limit of zero friction a — > 0, this action is reduced to the usual action for a charged particle in a 
homogeneous magnetic field 0. 

The set is linear one with constant coefficients. For such a case, the corresponding quantum Hamiltonian 
can be constructed according to ([62]) as 



Pi + Pi + 77 + QiPi + P2Q2 + ?2P 2 ) + PihQi ~ P1Q2) 



with operators qi and obeying canonical commutations relations, 

[qupj] = iSij, [qi,qj] = \pi,Pj] = 0, i,j = 1,2 
By help of a canonical transformation (pi, qi;p2, 92) (p> 2; 9> 2/)> where 



^ /i2 1 ~2\ 



(86) 



P = Pi + 7791; x = 9l> 9 = £»2 + 2 92, 2/ = 92 



we reduce (1861) to the form 



H = - 
2 



p 2 + g 2 + /3(gi-py) + 



(3 2 -a 2 



(87) 



Condition l|82[) in the case under consideration implies a <C /?, that is why a 2 will be omitted in (|87|l in what 
follows. 

Consider eigenstates \1/ for two mutually commuting operators H and L = py — qx, 

H^ = E^ , L<3> = . 
It is convenient to perform the following canonical transformation (fi, x; q, y) — > (P, X; Q, V), 



1 



P- -KV, X= - [q+-x 



/3 

<3 = *-§*,* = ± 



/3- 



0, 
P+ 2 y 



It is easy to see that 



H 



P 2 + (3 2 X 2 ) , L = P^iHi -H), Hi 



2V>2 



+ 0*Y 



Operators H and H\ are Hamiltonians of two independent harmonic oscillators. Then we can divide variables 
solving ([55)) . Thus, we obtain solution of the eigenvalue problem (|55|) . 

^ = ^ n , l {X,Y)=i> n {X)^ l {Y), E n = p(n + ~^ , M nl = l-n, n, I = 0, 1, 2, ... 

where ^ n and ^ are eigenstates of the Hamiltonians P and respectively (given e.g. by (|76[) ). Finally, 
stationary states ^(i) of the corresponding Schrodinger equation with the Hamiltonian H have the form 



-iE n t 



(89) 



We dehne the classical energy E of the system under consideration according to [33] as the mechanical 
energy of the system without friction, 



E 



p 2 +q 2 + (3(qx -py) + !L (x 2 + y 2 ) 



One can see that the energy depends of time as follows: E = E$e 2at . An operator E that corresponds to the 
classical quantity E reads: 



E 



1 



P l + f3 2 X 2 + a[XY - PQ 



(XY - PQ) 



-2at 



2a 



PY + QX 



e~ 2at + o 



12 



Mean values of this operator in stationary states (|89")1 can be easily calculated, they are 

Similar to the damped oscillator case considered above, at each fixed time instant, the energy spectrum is 
discrete, however, it decreses with time exactly as in classical theory. 

We would like to note that in the work fTU] it was shown that although an action principle for the second- 
order equations (|53")1 describing a radiating point-like charge does exist, none of the possible corresponding 
Lagrangians in the limit of a — > reduces to the Lagrangian of a particle in a magnetic field modulo a total 
time derivative. That is, in the case of a radiating point-like charge a perturbation (in the friction parameter 
a) of a second-order action does not correspond to a perturbation of the equations of motion ([53")) . For this 
reason, we expect some difficulties with the limit of a — > in the quantum theory of a radiating point-like 
charge resulting from quantization based on an action functional in the second-order form (such quantization 
for a damped harmonic oscillator was presented in |38j-[37j). 

8 Concluding remarks 

We stress that any nondegenerate set of differential equations written in an equivalent first-order form can 
be derived from an action principle. In the general case, such a set does not provide enough information to 
fix a class of quantum theories that, in the classical limit, provide this set of differential equations for mean 
values. Therefore, physical considerations must be used to choose an adeqate quantum theory. In particular, if 
one definitely knows that a non-Lagrangian set of equations describes a dissipative system, which is subjected 
to a dissipation due to essential interaction with an environment (reservoir), it is reasonable to consider the 
system and the reservoir as two interacting subsystems of a closed system. Then a quantum description of the 
dissipative subsystem can be obtained from a quantum theory of a whole system by averaging over the reservoir. 
Such an approach was developed in many articles, see |24j-|34j. However, one cannot consider such an approach 
as quantization of initial dissipative subsystem, since quantization was already made for the whole system. In 
the present article, we consider an approach where we actually quantize a system with a given set of equations. 
It turns out that its "non-Lagrangian" behavior is due to a time-dependent external field. It is a principally 
different physical situation in comparison with dissipation of a subsystem. However, quantum theories obtained 
from our procedure may be useful to describe some quantum-mechanical properties of both dissipative systems 
and "non-Lagrangian" systems of other physical nature, like a monopole. 

Acknowledgement 1 Gitman is grateful to the Brazilian foundations FAPESP and CNPq for permanent 
support; Kupriyanov thanks FAPESP for support. 

References 

[1] V.G. Kupriyanov, S.L. Lyakhovich and A. A. Sharapov, J.Phys.A38 (2005) 8039. 

[2] P.A.M. Dirac, Proc.Roy.Soc. A133 (1931) 60. 

[3] J. Douglas, Trans.Am.Math. Soc. 50, N 71 (1941) 71-87. 

[4] V.V. Dodonov, V.I. Man'ko and V.D. Skarzhinsky, Arbitrariness in the choice of action and quantization 
of the given classical equations of motion, preprint of P.N. Lebedev Physical Institute (1978). 

[5] W. Sarlet, J.Phys.A:Math.Gen.l5 (1982) 1503. 

[6] P. Havas, Actra.Phys.Aust. 38 (1973) 145. 

[7] R. Santilli, Ann.Phys. (NY)103 (1977) p.354-359. 

[8] S. Hojman and L. Urrutia, J.Math.Phys. 22 (1981) 1896. 

[9] M. Henneaux, Ann.Phys, 140 (1982) 45. 



13 



[10] V.G. Kupriyanov, Int.J.Theor.Phys.45 (2006) 1129. 

[11] M.Henncaux, L.Sheplley, J.Math.Phys 23 (1982) 2101. 

[12] J.Cislo and J.Lopuzanski, J.Math.Phys. 42 (2001) 5163. 

[13] P.Tempesta, E. Alfinito, R.Lco and G.Soliani, J.Math.Phys. 43 (2002) 3583. 

[14] D.M. Gitman and I.V. Tyutin, Quantization of Fields with Constraints (Springer- Verlag, Berlin 1990). 

[15] D.M. Gitman, S.L. Lyakhovich, M.D. Noskov and I.V. Tyutin, Izv VUZov 3 (1986). 

[16] D.M. Gitman and I.V.Tyutin, Nucl.Phys. B630 (2002) 509. 

[17] F.A. Bcrczin and M.A. Shubin, Schrodinger Equation (Kluwcr Publ., New York, Amsterdam 1991) 

[18] H.R. Lewis and W.B. Riesenfeld, J.Math.Phys. 10 (1969) 1458. 

[19] J.K. Kim and S.P. Kim, J.Phys.A32 (1999) 2711. 

[20] H.Kim and J.Ycc, Phys.Rcv. A66 (2002) 032117. 

[21] Kyu Hwang Yeon, Chung In UM, T.F.George, Phys.Rcv. A68 (2003) 052108. 

[22] H. Bateman, Phys.Rev. 38 (1931) 815. 

[23] E. Kanai, Progr.Theor.Phys., 3 (1948) 440. 

[24] I.R. Scnitzky, Phys.Rev. 119 (1960) 670. 

[25] R.P. Feynman and F.L. Vernon, Ann.Phys. 24 (1963) 118. 

[26] H.Haken, Rev.Mod.Phys. 47 (1975) 68. 

[27] H.Dckkcr, Phys.Rcv. A16 (1977) 2116. 

[28] A. Caldeira and A. Leggett, Physica A121 (1983) 587. 

[29] I.A.Pedrosa and B. Baseia, Phys.Rcv. D30 (1984) 765. 

[30] D.Walls, G.Miburn, Phys.Rev. A31 (1985) 2403. 

[31] H.Grabert, P.Shramm and G. Ingold, Phys.Rep 168 (1988) 115. 

[32] Li Hua Yu, Chang-Pu Sun, Phys.Rev. A49 (1994) 592. 

[33] I.Joichi, S.Matsumoto and MYoshimura, Prog.Theor.Phys., 98 (1997) 9. 

[34] M. Rosenau, A. Caldeira, S. Dutra, H. Westfahl, Phys.Rcv. A61 (2000) 022107. 

[35] W.E. Brittin, Phys. Rev. 77 (1950) 396; N.A. Lcmos, Phys.Rev. D24 (1981) 2338. 

[36] R.Hasse, J.Math.Phys. 16 (1975) 2005. 

[37] D.C.Khandckar and S.V.Lavandc, J.Math.Phys. 20 (1979) 1870. 

[38] H.Dekker, Phys.Rep. 80 (1981) 1. 

[39] E. Celeghini, M. Rasetti and G. Vitcllo, Ann.Phys. 215 (1992) 156. 

[40] R.Banerjee and P.Mukhcrjee, J.Phys.A35 (2002) 5591. 

[41] M. Blasonc and P. Jizba, Ann.Phys. 312 (2004) 354. 

[42] D. Chruscinski, Ann. Phys. 321 (2006) 840; Ann. Phys. 321 (2006) 854 

[43] L.D. Landau and E.M. LIfshitz, Mechanics, Oxford:Pcrgamon, 1962. 

[44] V.L. Ginzburg, Theoretical Physics and Astrophysics: Aditional Chapters, Moscow: Nauka, 1987 (in 
russian). 



14 



